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Abstract 

The two-loop ultrasoft contributions to the next-to-leading logarithmic (NLL) running of the 
QCD potentials at order v'^ are determined. The results represent an important step towards 
the next-to-next-to-leading logarithmic (NNLL) description of heavy quark pair production and 
annihilation close to threshold. 
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I. INTRODUCTION 



The measurement of the hne-shape of the total top-antitop quark production cross section 
in the threshold region ^/s ~ 2mt is one of the major tasks within the top quark physics 
program at the ILC. The most prominent quantity to be measured is the top quark mass, 
and one expects an improvement in precision of rrit by about an order of magnitude to mass 
measurements based on reconstruction obtained at the Tevatron and the LHC [1]. 

To obtain a meaningful theoretical description of the nonrelativistic threshold dynamics 
it is required to systematically sum the so-called Coulomb singular terms oc (as/f)" in a 
systematic nonrelativistic expansion, v being the relative velocity of the top quarks. This 
task is achieved by means of effective theories based on nonrelativistic QCD (NRQCD) 0]. 
In this approach, however, sizeable logarithmic terms oc (a^lnf)" are not systematically 
accounted for, which leads to rather large normalization uncertainties of the cross section 
line-shape. The next-to- next-to-leading order predictions in this "fixed-order" approach 
were estimated to have a normalization uncertainty of order 20% p|. These normalization 
uncertainties do not affect the top mass measurement, which primarily depends on the cm. 
energy where the cross section rises, but they render measurements of other quantities such 
as the top total width or the top quark couplings impossible. To match the statistical 
uncertainties, that are expected for these quantities, a theoretical precision of the cross 
section normalization of at least 3% would be required. 

In Refs. [^] (see also Ref. [sj) it was demonstrated that the summation of logarith- 
mic (aslnf)" terms, using renormalization-group-improved perturbation theory, can sig- 
nificantly reduce the normalization uncertainties of the threshold cross section. Concerning 
QCD effects the renormalization group improved LL (leading-logarithmic) and NLL order 
predictions of the threshold cross section are completely known, but no full NNLL order 
prediction exists at present. The full NNLL order prediction is, however, required to obtain 
a reliable estimate of the remaining theoretical normalization uncertainties. 

The missing ingredient is the NNLL running of the Wilson coefficient of the leading order 
effective current, that describes production and annihilation of a nonrelativistic tt pair in a 
S-wave spin-triplet state ('^5'i). Adopting the label notation from vNRQCD [6, 7] the current 
has the form 



Ji,p = i'lcr{ia2)x* 



p ' 



(1) 



where ipp and Xp annihilate top and antitop quarks with three-momentum p, respectively, 
and where color indices have been suppressed. The current does not have a LL anomalous 
dimension because there is no one-loop vertex diagram in the effective theory that contains 
UV divergences associated with the current Ji^p. Such UV divergences arise at NLL order 
from insertions of the NNLL kinetic energy operators and from insertions of the NNLL order 
potentials |6|. The corresponding computations were carried out in Refs. 0, B [o], [13] and are 
completed. Using the conventions from j3| the resulting NLL order renormalization group 
equation for the Wilson coefficient ci of the current has the form (S^ = 2) 



d 

'^^ln[ci(z/)] 



-I- a^^mu) 



V) 



^{Cp-2Ca) 



V] 



(2) 
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where u is the vNRQCD velocity renormahzation parameter that is conveniently used to 
parametrize the correlation between soft and ultrasoft dynamical scales within the renor- 
malized effective theory [6|. The analogous evolution equation for currents describing pairs 
of quarks and colored scalars in any angular momentum and spin state (^'^^^Lj) were de- 
rived recently in Ref. 11 1. In Eq. ([2]) the term v'f^ is the Wilson coefficient of the Coulomb 
potential oc 1/k^, and V2''^ and Vr^^ are the coefficients of the 0{v'^) potentials with the 
momentum structure 1/rn? and (p^ + p'^)/(2m^k^), respectively, m being the heavy quark 
mass. The coefficients V^^^^g from the so-called non-Abelian potentials that scale like 

l/(m|k|), and vi^'^ is the coefficient of the spin-dependent potential that can contribute for 
spin triplet S-wave states. The superscripts (s) refer to the color singlet state of the quark 
pair. 

At NNLL order there are two kinds of contributions to the evolution of Ci that have 
to be accounted for. The first kind arises from three-loop vertex diagrams that come from 
insertions of subleading soft matrix element corrections to the potentials and from insertions 
of potentials with additional exchange of ultrasoft gluons. The corresponding computations 
were carried out in Ref. 13] ■ These effects are referred to as the non-mixing contributions 



as they affect the evolution of ci directly. The second kind of contributions arises from 
the subleading evolution of the potential Wilson coefficients that appear in the NLL or- 
der renormahzation group equation shown in Eq. ([2]). They are referred to as the mixing 
contributions as they affect the evolution of Ci indirectly. Except for the coefficient of the 
Coulomb potential V'f'^ 0, [III and for the spin-dependent potential vi'^'' [3] no complete 
determination for the subleading evolution exists at present. The analysis of the three- loop 
(non-mixing) terms in Ref. [13] revealed that the contributions involving the exchange of 
ultrasoft gluons are more than an order of magnitude larger than those arising from soft 
matrix element insertions and are as large as the known NLL contributions. The reason is 
related to the larger size of the ultrasoft coupling as{mv'^) and to a rather large coefficient 
multiplying the ultrasoft contributions. These large ultrasoft contributions are responsible 
for an uncertainty in the normalization of the most up-to-date threshold cross sectionpre- 
diction of at best 6% [l^, which is quite far from the required precision. (See also Ref. [sf for 
an alternative analysis without NNLL non-mixing effects.) The corresponding effects of the 
soft matrix element corrections are only at the level of several per mille. From this analysis 
it is reasonable to assume that the ultrasoft effects, which form a gauge- invariant subset, 
also dominate the mixing contributions. This assumption is also consistent with the results 
for the NLL evolution of the coefficient vi*^ [13], which is dominated by soft effects, since 
the ultrasoft gluon coupling to heavy quarks is spin-independent at this order, and which 



was found to have a very small numerical effect as well [16|. From the parametric point of 
view this can be understood from the fact that ultrasoft effects can affect the coefficients of 
the spin-dependent potentials only indirectly through mixing via potential loop divergences, 
since at this order the coupling of ultrasoft gluons to heavy quarks are spin- independent. For 
the coefficient of spin-dependent potentials this leads to ultrasoft NLL terms oc a|(a;7lnf )", 
where as and au are the strong coupling at the soft and ultrasoft scales, respectively. These 
terms are parametrically suppressed compared to ultrasoft NLL terms oc asaui^au^^vY 
that can affect the spin-independent potential coefficients through two-loop ultrasoft loop 
diagrams. We can therefore expect that ultrasoft effects are dominating the higher order 
evolution of the spin-independent potentials. 

In this work we determine as a first step the two-loop ultrasoft contributions to the NLL 
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anomalous dimensions of the spin- independent potentials V2 and Vr- The outline of this 
work is as follows: in Sec. [TTl we briefly review the effective theory setup used for our work 
and in Sec. [TTl]we present our results. Sec. [IV] contains a brief numerical analysis and Sec. |V] 
the conclusion. 



II. THEORETICAL SETUP 

The vNRQCD Lagrangian basically consists of three parts jl, 0], 

■^vNRQCD = + + , 



(3) 



containing kinetic terms and ultrasoft interactions, potential interactions and interactions 
involving soft degrees of freedom, respectively. The "ultrasoft part" has the form 



r 



p 



2m 



P 

Sm^ 



+ . . . 



-G^^G, , + (4) 



where the ultrasoft gauge- covariant derivative reads = d'^ + iguA^{x) and gu is the 
strong coupling at the ultrasoft scale mz/^. The terms describe besides the propagation of 
the heavy quarks their interaction with ultrasoft gluons. The "potential part" Cp describes 
the potential quark- ant iquark interactions and reads 



- ^"/^Ar (P, P) ^p/„ X-p'A X-pr + 

P.p' 

where a, /3, A, r are color indices and 



(5) 



+ 



+ (1 



«/3 



LArJ 



k2 



mlkl 



mlkl 



+ 



+ 



Vf^(p^ + P^^) 

Vr^'^(p^ + P^^) 



V. 



(T) 



+ 



+ 



V. 



(1) 



m^ 



+ 



(6) 



with k = p' — p. The spin-dependent potentials have not been written out in Eq. IQ, 

since they will not be relevant in this work. One can convert the potentials from the basis 
formed by the two 3 ® 3 color structures 1 (S> 1 and T"^ (8> into the more physical color 
singlet /octet basis by 



kinglet 




K)ctet 





1 



Vt(S)T 



(7) 



In the soft sector of the theory there are 6-field operators needed for the description of 
soft Compton scattering off quark- ant iquark potentials (see Fig. [1^). As was pointed out in 
Ref. [7] these operators have zero matching conditions at the hard scale u = 1 and contribute 
to the ultrasoft renormalization of the 0{v'^) potentials. The operators have the structure 



^2c 



9t 



(8) 
(9) 
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FIG. 1: a) Operator 02i, b) operator with soft lines closed to a loop, c) diagram with soft and 
ultrasoft loop contributing to potential counterterms. 



where (pg, and Cq denote soft quarks, gluons and ghosts respectively, gs is the strong 
coupling at the soft scale mz/ and a is the order at which the operators contribute in the 
f -counting. For the explicit form of the the F's we refer to Ref. [7(] (see also Ref. 13]). The 



superscript (T) refers to the color structure (T^ ® T^) that arises upon closing the two 
soft lines (see Fig. [T]d). There is another set of such operators, which are not displayed, 
having the superscript (1), that refers to the^ color structure (1 1) that arises upon closing 
the two soft lines. As was shown in Ref. [tJ, closing the two respective two soft lines for 
the leading a = operators, one obtains a structure for the one-loop four-quark matrix 
elements that is identical to the one-loop soft-matrix element corrections to the Coulomb 
potential. Moreover, the Wilson coefficients of the a = 2 operators become nonzero for 
V < 1 from ultrasoft UV divergences (and the associated pull-up terms 0, [l8|) in higher- 
order corrections to the vNRQCD matrix elements describing Compton scattering off a 
quark- ant iquark pair. The coefficients of these operators are directly affected by ultrasoft 
gluon corrections to this process. Through the UV divergences that arise upon closing the 
two soft lines they mix into the potential coefficients and Vt^\ Using the counting 
mentioned above the contributions from these mixing contributions to the evolution of V2*^ 

is] 

and Vr are parametrically suppressed compared to the contributions from the two-loop 
ultrasoft diagrams that renormalize these potentials directly. For the purpose of this work 
we need the a = 2 operators that have the form 

,^(2),(i) _ n^2),(T) _ n^o),{T) 

i i 

og"" = E OS"" . og'-'-' = E OS""' . (10) 

i i 
(\ T) (1 

which have the Wilson coefficients Cg;^' and C2p' , respectively. We note that the mixing 

(s) (s) 

effects of these operators lead to numerical effects for V2 and Vr that are substantially 

(s) (s) 

smaller than those from the two-loop ultrasoft diagrams that renormalize V2 and Vr 
directly. We include them because their contribution leads to a rather simple form for the 
two-loop ultrasoft corrections to V2'^'' and Vr*"* we determine in this work and also in order 
to illustrate the size of the contributions not considered in this work. 

All fields, couplings and Wilson coefficients in the above Lagrangian are to be understood 
as bare quantities. For the renormalized ones we chose the usual conventions in d = 4 — 2e 



dimensions {fis = mi^, fJ'U = Tnv"^): 

For convenience, we will drop the index R throughout this paper and only deal with MS 
renormalized quantities in the following. 



III. RESULTS 



In general the counter term of the potentials induced by ultrasoft renormalization takes 
the form 

6V = AV + C, (12) 
where 5V , V and C are 2-vectors and A is a 2 x 2 matrix in the (1 ® 1, T"^ ® T"^] 



V 



Vim 



space, 
(13) 



From Eq. f|T2l) one can deduce the general form of the renormalization group equations for 
the potentials: 



u-^V = -2tV -{1 + AY^{2tC + u-^C + {u-^A)V) . 

/~1 T I ^ /~1 T I /-ill ' 



. _ , (14) 

du dv dv 

The full one-loop results for A and C can be found in Ref. 0, [1]. For the ultrasoft 
contributions at the two-loop level we define different classes of (divergent) two-loop 
diagrams in Feynman gauge contributing to the matrix A: 



class topology insertions/vertices (e.g. ip'^pgA^ipp, tp^g^^tpp, ■ip'^^^ipp, etc.) order 



1 

2 
3 
4 
5 



8 
9 



Fig 
Fig 
Fig 
Fig 
Fig 



[2^ four gA^ vertices 

one gluon selfenergy insertion and two gA^ vertices 

[2^ two gA^ and two g^^ vertices 

[2}3 one gluon selfenergy insertion and two g^^ vertices 

[2H two gA'^ vertices and two insertions of the ^ operator on the 0{v'^) 
internal quark lines 

Fig. [2)d one gluon selfenergy insertion, two gA^ vertices and two inser- (9(t'^) 
tions of the operator on the internal quark lines 

Fig. [2)d one gluon selfenergy insertion, one gA^, one g^^ vertex and one 0{v'^) 
insertion of the operator on the internal quark lines 

Fig. ^ one triple gluon interaction, one gA^ and two g^^ vertices 0{v'^) 

Fig. [2t one triple gluon interaction, two gA^, one g^^ vertex and one 0{v'^) 
insertion of the operator on the internal quark lines 
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To each of these classes there is a corresponding class of diagrams, which contribute to the 









FIG. 2: 2-loop topologies: a) "Abelian", b) with gluon selfenergy, c) "Non-Abelian" . Left-right 
and up-down mirror graphs are not shown. 



heavy quark wavefunction counterterm SZ^^p, with the same type and number of vertices 
and insertions. The results for the different contributions to SZ^^p after subtraction of the 
respective one-loop sub divergencies are shown in Tab. [H The table also shows an example 
diagram for each of the classes. The result for the ultrasoft two-loop contributions to the 
matrix A reads: 



A = ^ 



_|_ "(7 Ca 



2,2 



k2 (C^+^_3^)k2 + C^(p2 + p'2) 



Ci(k2-f(p^ + p'^)) 



p 



(15) 



where Ci = - 2Cf) and Cd = 8Cf - SCa- 



Here and in the following au = a<j(mz/^) and as = ^^(mz/). The expression in Eq. f|T5|) 
represents the main result of this work. It is an important cross check that the (9(a^f°) 
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class 


example diagram 


contribution to 6Z^^p 


order 


1 




Cf{Cf-Ca)»1 I CACf-af^ 




2 




CF(5CA-4Tn/)a2^ CF{13CA-8Tn^)a2^ 
327r^€^ ' 487r^e 




3 












4 





CF(15CA-12Tn/)p2Q:2^ Cf (46 Ca -32 Tn/)p2a2^ 
288 m'^TT-^e-^ 288m'^7r-^e 




5 


-^-® — ^— ^ — 


CFiCF-CA)p''al CACFY>^al 




6 




CF(75CA-60Tn/)p2a2^ CF(110CA-64rn/)p2a2^ 




— ^^^s g^^^ 


288m27r2e2 ' 288™^^^^ 


7 




Cf (15 CA-12Tn/)p2a2^ Cf (16 Ca -8 Tnj)p2a2^ 






72m^7r2e^ 72m^7r2e 
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TABLE I: Contributions to 6Z^,^p at ultrasoft 2-loop level {ajj = as{mv'^)). 



terms vanish identically as a consequence of ultrasoft gauge invariance. In addition, our 
result for the ultrasoft two-loop 0(1;°) wave function renormalization constant agrees with 
the 0{a'^m^) heavy quark wave function renormalization constant in heavy quark effective 
theory (HQET), see e.g. Ref. [l9|. There is agreement since the purely ultrasoft sector 
of NRQCD is closely related to HQET [20|. Furthermore we find that all (^(f^) terms 
cancel in our result for the ultrasoft two-loop wave function renormalization constant (see 
Tab. [I]). This cancellation also takes place for the 0(f 2) terms at one-loop and is required by 
consistency since the form of the bilinear heavy quark kinetic terms in Eq. (jlj) is protected 
from quantum corrections by reparametrization invariance. As an additional cross check we 
have also determined all two-loop fermionic Tuf corrections in Coulomb gauge [2l|. 

From the results obtained in Eq. f|T5|) one can deduce without effort the ultrasoft con- 



(1 T) 

tributions to the NLL evolution of the Wilson coefficients C2i^2pt since the additional soft 
external lines describing soft Compton scattering do not affect the structure of the ultra- 
soft loop diagrams. Including the two-loop soft terms induced by the pull-up mechanism in 
analogy to the 1-loop case [?, 1^ and the corresponding (known) one- loop contributions the 
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respective renormalization group equations read 

2 4^'^) (2a^-a5)+4i?(i'^) {2al - 



d ^{i,T)/ X 



d ^(i,T)/ N 



a 



s) ? 



where the constants A\^'^ and ^j^p are defined as 



(16) 
(17) 



4) = 0. 



16 



AT 



A 



(1) 



A 



(T) 



Ca 

67r 
2^ ' 



p 

B^P = Ca 



CaCf{Ca — 2Cf) 



C.(-577r2 + 188) - 40Tnf SC.C^ 



5, 



(1) 



4327r2 



Ci?(C^ — 2Cp 



- H 

12 



Ca{21tx^ - 47) + l^Trij 
' 2167r2 



(T) 



(Ca — 2Cf 



CA(37r^ -47) + lOTnj 
' 72^^2 



The solutions respecting the matching conditions C^2fc 2p = at z/ = 1 have the form 



(/3o = ¥C.4 - iTn, , ft = f d - 4CFrny - f C^Tn^) 



AirA 



AttA 



au 



— In 



Ai'^^f3, 




PI 


Po 




8vr4^'^) 


PI 


Po 



[au - as) 
{au - as) 



(18) 



(19) 



Upon closing up the two soft gluon, ghost and hght quark hues one arrives at the following 
mixing contributions for C in Eq. f|T2l) : 



C 



asPo 
' k2e 



^(1) ki , ^(1) p^+p'" 



It is then straightforward to derive the evolution equation for the coefficients and V, 
of the spin-independent (9(f^) potentials. Using 



(20) 



V 





k2 



(21) 



with Vc(z^) = 4 TT asimu) on the RHS of Eq. ( |T^ . and including again the soft pull-up terms 
we obtain 



— V- 

d_ 
dv 



87r A^''^) as (2ac/ - a^) + 167r ^^'^^ (2a^ - a|) - 2 /3o a| C^^'^^ 



V(^''^) (z/) = 87r «5 (2«c/ - «5) + 16vr M^'^) (2«^ - «|) - 2 /3o «| C^^''^) , (22) 
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for the dominant NLL ultrasoft contributions. To render the anomalous dimensions finite, 
we also have to account for UV-divergent contributions from diagrams of the type in Fig. [1]:. 
Note that in Eqs. fl22|) we have also displayed the LL ultrasoft terms including the respective 
soft pull-up terms. The full set of the known LL contributions from soft one-loop diagrams 
has not been displayed. 

From the previous equations we find the following simple result for the ultrasoft one-loop 
and two-loop contributions to the evolution of the coefficients V2 and Vr up to NLL order: 



Ait aJmu) C, 



V] 



(23) 



The corresponding coefficients for the color singlet potentials read 

V^'^M)^^ = 47ra,(mz.) (c«(z/) - C^Cf ^(z/)) , 

'yl'K^))^^ = Snasimu) (c«(z.) - C^Cf )(z.)) . (24) 

It is now straightforward to derive from Eq. ([2]) the two-loop ultrasoft part of the NNLL 
mixing contributions to the running of the ^Si current coefficient Ci from the coefficients 
V2'^^ and vi^'' . Parametrizing the form of ci as [iT 



In 



ci(l) 



(25) 



where i^^^^^ and ^nm^^ refer to the NNLL mixing and non-mixing contributions, respectively, 
we find 



NNLL 
m;2r,usoft 



^^^"sM 

Po 



7 TT^ 



z 1 - In 



2-z 



\nz 

4 2 



-In^ 



2-z 



2Li9 



Stt'Cf 5 2 



Po 



3-2z- z^ -A ln(2 - z) 



(26) 



where 



A = A«+2A«-C^(Af + 24^) 
B = i?(^)+2i?«-C^(5f) + 2i?f)), 



Oisirn) 



LL 



1 + ^i^^ln. 
27r 



(27) 



For completeness we also give the form of the ultrasoft corrections to from V2 and Vr-: 



fNLL ^ StT^C p ~ 
?2r,usoft n2 

Po 



Aasim) -1 + z + {2- z)ln{2- z) 



(28) 
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IV. NUMERICAL DISCUSSION 



In Fig. [3^,b the renormalization parameter v dependence of the color singlet coefficients 
and Vr*"* in Eg. (12^ is shown at LL order (dashed lines) including all soft and ultrasoft 
contributions 0, [13] , and in addition including the full two-loop ultrasoft contributions with 
(solid lines) and without the mixing contributions from the zero matching operators (dotted 
lines), that are represented by the respective last terms on the RHS of Eqs. fl2^ . For the 
heavy quark mass we have adopted rrit = 175 GeV for the application to the top-antitop 
quark threshold, and as{mt) = 0.1074. To further illustrate the importance of the two-loop 




FIG. 3: Renormalization scale dependence of V2 and Vr ■ In a), b) the full LL results (soft 
+ ultrasoft) (dashed lines), plus the NLL ultrasoft contributions with (solid lines) and without 
(dotted lines) the mixing terms is shown. In c), d) only the corresponding ultrasoft results are 
displayed omitting the known soft LL running. All of them include soft pull-up terms. 



ultrasoft contributions we have displayed in Fig. [3b, d the evolution of the coefficients at 
LL order due to the one-loop ultrasoft contributions alone (dashed lines), and in addition 
including the two-loop ultrasoft contributions analogous to Figs. [3^ and b. Comparing the 
behavior of the LL order and NLL two-loop ultrasoft contributions in Fig. we find 
that the NLL order corrections are of the same size as the ultrasoft LL order contributions. 
While for the potential this is of no concern, since for V2''^ the ultrasoft contributions 
are suppressed by a small overall color factor compared to the soft LL contributions, the 
result for Vr^^ shows that the two-loop ultrasoft corrections are indeed anomalously large, 
similar to the NNLL ultrasoft non-mixing contributions to the evolution of Ci determined 
in Ref . T3] . From the dotted lines in Fig. [3] we also see that the correction arising from the 
mixing contributions are numerically small, which is consistent with the parametric counting 
mentioned in the introduction. 
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It is also instructive to analyse the impact of the NLL two-loop ultrasoft contributions in 
and Vr*'* on the NNLL evolution of the ^Si current coefficient ci. Evaluating Eq. (!26|) 
we find {6ci)^^^^''^'' ^ ^m>,usoft = (-1-9%, -0.5%) for u = 0.1,0.2. These corrections lead 
to a partial compensation of the anomalously large ultrasoft corrections that dominate the 
NNLL non-mixing contributions to the evolution of Ci. It remains to be seen whether the 
two-loop ultrasoft corrections to the evolution of the l/(m|k|) QCD potentials will lead to 
a similar effect. We will address this issue in a subsequent work. 



V. CONCLUSION 

We have determined the two-loop ultrasoft contributions to the NLL order renormal- 
ization group equations of the spin-independent QCD potentials V2 (1/m^) and Vr 

((P^ + p'^)/(2m^k^)). We find that the two-loop ultrasoft corrections are larger than the 
ultrasoft contributions at the one-loop level and lead to a partial cancellation of the anoma- 
lously large NNLL ultrasoft non-mixing contributions to the renormalization group evolution 
of the leading ^5*1 current that describes top-antitop pair production at threshold in e+e" 
annihilation. 
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